Transport properties of molecular junctions from many-body perturbation theory 
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The conductance of single molecule junctions is calculated using a Landauer approach combined 
to many-body perturbation theory (MBPT) to account for electron correlation. The mere correc- 
tion of the density-functional theory eigenvalues, which is the standard procedure for quasiparticle 
calculations within MBPT, is found not to affect noticeably the zero-bias conductance. To reduce 
it and so improve the agreement with the experiments, the wavefunctions also need to be updated 
by including the non-diagonal elements of the self-energy operator. 

PACS numbers: 85.65.-|~h,71.10.-w,72.10.-d,73.40.-c,73.63.-b 



Most recent theoretical studies of coherent transport 
in nanojuctions are based on a Landauer approach [l| , in 
which the electron interactions are treated at a simplified 
mean- field level using density- functional theory (DFT). 
While this approach has proven quite successful for sys- 
tems having a strong coupling between the molecule and 
the metallic leads 0, it overestimates the zero-bias 
conductance of weakly coupled systems (by up to 3 or- 
ders of magnitude) compared to experimental measure- 
ments l^-la]- Many explanations have been proposed for 
such a discrepancy. For example, arguing an uncertainty 
over the experimental junction structure, the sensitive- 
ness to the contact geometry was investigated 0, Q ■ 

More fundamentally, it has been found that the zero- 
bias conductances may vary by several orders of mag- 
nitude when using different exchange-correlation (XC) 
energy functionals i, 113 . Moreover, the DFT lack of 
derivative discontinuity was shown to be a source of sig- 
nificant errors in weakly coupled systems In fact, the 
validity of using DFT to treat electron interactions in the 
Landauer formalism has also been questioned 



it is believed that a many-body theory should be better 
suited. Indeed, the DFT eigenstates do not have a formal 
direct connection to the quasiparticle (QP) states. 

Recently, many-body perturbation theory (MBPT) 
calculations [l3j have been used in combination with the 
Landauer formalism to study the conductance of a sim- 
ple gold monoatomic chain [ij] . The applicability of this 
method to more complex systems, as realistic three di- 
mensional junctions, is still out of reach. Therefore, a 
simpler approach relying on a model self-energy [l^ , 
based on projectors onto molecular orbitals, has been ap- 
plied to molecular junctions with some success. Never- 
theless, a formal justification of this model is still lacking. 

In this letter, the effect of electron-electron interac- 
tions on the zero-bias conductance is investigated for 



the benzene diamine (BDA) and benzene dithiol (BDT) 
molecules attached to gold electrodes. The true QP elec- 
tronic structure, calculated using different approxima- 
tions for the MBPT self-energy, is used in the Landauer 
formalism. It is found that correcting only the eigenener- 
gies has no or little effect on the conductance. To improve 
the agreement with the experiment, it is thus crucial to 
update also the wavefunctions by taking into account the 
off-diagonal self-energy matrix elements. A comparison 
of the original and the updated local DOS (LDOS) at the 
Fermi level shows that both the molecule and the gold 
regions are affected. In particular, it is found that the 
change of the zero-bias conductance is triggered by a de- 
crease of the molecular character and an increase of the 
eg{dz2) character of the wavefunctions on the gold atoms. 
Finally, our results arc compared to those provided by the 
molecular-projectors model proposed in Refs. 
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The DFT and MBPT calculations arc performed using 
the Abinit package [l^ . The XC energy is approximated 
by the PBE functional [l3| • Norm-conserving pseudopo- 



1^, and tentials 18[ are used. For gold, these include the 5s 



and 5p semicore states which are crucial for the MBPT 
calculations [l9j . The wavefunctions are expanded on 
a plane-wave basis set up to kinetic energy cutoff of 
30 Ha. The corresponding maximally-localized wannier 
functions (MLWFs) [2^ are obtained following the proce- 
dure of Rcf. [2l| . These are used in the transport calcula- 
tions which are performed with the WanT package p^ . 

Two systems are considered: BDA and BDT attached 
to gold electrodes |23| . For the atomic relaxation, a 2 x 2 
surface cell is adopted for the Au (111) surface with seven 
atomic layers in the electrodes. The relaxed geometries 
agree with previous theoretical predictions [s, 15, 24|. For 
the MBPT calculations, three gold layers are removed. 
This guarantees an affordable computational cost with- 
out compromising the calculated conductance [25| de- 
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spite the rather packed geometry (with ~3.1 A between 
the repeated images of the molecules). All the param- 
eters of the calculations are converged to ensure an er- 
ror smaller than 0.001 Qq on the zero-bias conductance 
G{E=0). Hence, a 8x8x3 grid of k points is adopted to 
sample the Brillouin zone. The QP corrections are calcu- 
lated explicitly for 210 bands at 96 irreducible k-points 
including ^^300 bands in the calculations. 

For the MBPT calculations, the GW approximation is 
adopted for the self-energy S, neglecting vertex correc- 
tions 13[. In principle, E is a non-hermitean, non-local 
and frequency-dependent operator. The non-hermitian 
part of S is ignored, assuming infinite QP lifetimes; while 
the non- locality is fully taken into account. Different fla- 
vors arc then considered for the GW self-energy. In the 
standard single-shot GqWq approach, the self-energy is 
approximated using the DFT electronic structure; and, 
the QP corrections to the DFT eigenvalues are deter- 
mined using first-order perturbation theory considering 
only the diagonal elements of S in the space of the DFT 
orbitals (pl^^ of the contacted-molecule system. The fre- 
quency dependence is obtained using a plasmon pole 
model [2^ . In order to calculate the first order correction 
to the DFT wavefunctions, the off-diagonal elements of 
S would also be needed. Currently, such calculations are 
out of reach for the systems considered here. Neverthe- 
less, a full diagonalization can be performed if a further 
simplified self-energy is considered. For this purpose, the 
Coulomb- hole screened-exchange approach 27[ , which is 
a static approximation to the GW self-energy, is adopted 
here. For the sake of comparison, both the full (CHSX) 
and diagonal (CHSX) forms of this self-energy are con- 
sidered. Finally, the molecular-projectors model (MPM) 
proposed in Rcfs. [H, [l^ is also used for the self-energy: 
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where m runs over the DFT orbitals (/)™°' of the isolated 
molecule, Ep is the Fermi energy, and A is the scissor- 
operator shift to be applied. Its value is obtained in a 
two-step procedure. First, the correction to the DFT 
HOMO-LUMO gap of the molecule in the gas phase is 
computed either from a ASCF calculation or, more sim- 
ply, by matching the experimental value. Second, the 
presence of the metallic surface is taken into account 
through a classical image-charge model leading to a re- 
duction of A compared to the gas phase, as calculated 
using the GoWo approximation 28|, l29|. This self-energy 
operator is non-diagonal in the space of the DFT orbitals 
(p^^^ of the contacted-molecule system. For the sake of 
comparison, its diagonal form (MPM) is also considered. 

The values of the zero-bias conductance Q{E=0) cal- 
culated with these different approaches are reported in 
Table U for BDA and BDT. For the latter, the complete 
energy dependence of the conductance GiE) is also given 
in Fig. [TJ The DFT results are in good agreement with 



previous calculations (i-H, i, [l!] . Quite surprisingly, the 
Go Wo results for Q(E=0) are almost identical to the DFT 
ones despite the changes in the curve G{E). 
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TABLE I. Zero-bias conductance G{E=0) in Go units calcu- 
lated with diflferent approaches (see text) for the two systems. 

In order to understand this finding, the projected den- 
sity of states (PDOS) on the molecule is computed as 
reported in the bottom panel of Fig. [l]for BDT. In par- 
ticular, an inset shows a decomposition of the total PDOS 
on the molecule calculated within DFT in terms of the 
contributions of the different molecular orbitals close to 
the Fermi level. The HOMO-1 and LUMO (in yellow 
and pink, respectively) lead to very sharp peaks in the 
PDOS. In fact, these states are very localized on the 
molecule and present very little hybridization with the 
leads. Hence, they do not play an important role on the 
conductance at eV. The latter is thus mainly driven by 
the HOMO and LUMO-l-1 (in orange and cyan, respec- 
tively) in agreement with the findings of Ref. 24 1. 

As expected, the GoWo corrections open the gap be- 
tween the occupied and unoccupied molecular-like levels, 
as can be seen in the PDOS (bottom panel) and the con- 
ductance (top panel) in Fig. [T] This opening of the gap 
is much smaller than the one obtained for the isolated 
molecule (~5.5 eV). In particular, the gap between the 
HOMO and LUMO-Hl is increased only by ~0.5 eV. The 
presence of the metallic surface can be estimated to ac- 
count for a reduction of ^2 eV with respect to the iso- 
lated molecule H, 2^. More importantly, the GqWq cor- 
rections depend on the weights of the different states of 
the system on the molecular orbitals wt = \{(t>T^\(t^m°^)\'^ ■ 
These are quite small for the very hybridized HOMO 
and LUMO-l-1, hence the small gap opening. In con- 
trast, these weights are higher for the more localized 
LUMO and HOMO-1. Thus, the correction to their gap 
('--'1.2 eV) is more than twice bigger than for the HOMO 
and LUMO-l-1. Similar results are observed for BDA. 

Due to this very limited opening of the gap, the result- 
ing Go Wo zero-bias conductance Q{E=0) is almost equal 
to DFT one for both systems. A similar result is observed 
when using the CHSX approximation, even though this 
approximation usually leads to bigger gaps than those 
obtained with the Go Wo approach. In conclusion, the 
conductance G{E=0) is almost not affected compared to 
DFT when the self-energy only corrects the eigenvalues. 

When updating not only the eigenenergies but also the 
wavefunctions by the CHSX approximation, Q{E=0) is 
reduced by a factor 1.4 (1.7) in the direction of the exper- 
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E - Ep [eV] 

FIG. 1. Conductance (top) and PDOS on the molecule (bot- 
tom) as a function of the energy (the Fermi energy Ef is set 
to zero) for the BDT junction. Different approaches are con- 
sidered: DFT (in bl ack), G qWq (in blue), CHSX (in purple), 
CHSX (in red), and MPM (in green). For each panel, an in- 
set provides a zoom around the Fermi energy. In the bottom 
panel, a second inset shows the total PDOS on the molecule 
computed in DFT (black line), and the corresponding contri- 
butions of the HOMO-l (in yellow), the HOMO (in orange), 
the LUMO (in pink), and the LUMO-l-1 (in cyan). 

imcntal value reaching 0.013 (0.020) for BDA (BDT). 
The agreement with experiments would probably be fur- 
ther improved if one could afford to perform calculations 
with (i) more than one step in the self-consistency loop 
on the wavefunctions; (ii) more realistic geometries, such 
as those depicted in Ref. 0| ; and (iii) a higher number of 
gold layers entering the Landauer formula. 

The PDOS on the molecule (Fig. [J) shows that CHSX 
corrections contribute to further separate the unoccupied 
and occupied molecular-like orbitals. Since the eigenen- 
ergies are only affected at the second order compared 
to CHSX, this increased opening of the gap is to be 
attributed to modifications of the metal-molecule hy- 
bridizations. To gain more insight on the changes of the 
wavefunctions, the LDOS is computed in an energy win- 




FIG. 2. Difference between the CHSX and CHSX LDOS cal- 
culated in an energy window of 0.8 eV around the Fermi level 
for BDA (top) and BDT (bottom). Four isovalues are repre- 
sented: +4p in orange, -f-lp in dark red, — Ip in dark blue, 
and -4p in light blue, with p = 4 x 10"* (6 x 10"") e'/A^ in 
BDA (BDT). The Au, C, N, S, and H atoms are represented 
by yellow, gray, brown, green, and white spheres, respectively. 

dow of 0.8 eV around the Fermi level Ep'- 
LDOS(r)=/ y2\cl)f%r)\H{E~E,)dE. (2) 

JEp-OAeV J 

The effect of the off-diagonal elements of E can be ana- 
lyzed by plotting the difference between the CHSX and 
CHSX LDOS, as shown in Fig. [5] Two main changes can 
be identified. First, the molecular character is reduced. 
This charge transfer does not take place in space (it is less 
0.02e" from the molecule to the gold leads) but in en- 
ergy: the molecular orbitals being shifted away from Ep 
(see the PDOS in Fig. [T]). Second, the eg((ij,2) character 
increases on gold atoms (see the typical ring shaped red 
lobes on the gold atoms). For the Au add-atoms, the z 
direction of the eg{dz2 ) orbital is oriented along the Au-N 
bonding direction. For the atoms of the next gold layer, 
it changes though it is still conditioned by the Au-N di- 
rection. Since there is not a noticeable reduction of the 
gold character in other directions and no molecule-gold 
charge transfer, this increase of the eg{d^2) character can 
be mainly to a transfer of electrons from lower energies to 
the Ep region due to the complete diagonalization of the 
self-energy corrected hamiltonian. In contrast, the dif- 
ference between the DFT and CHSX LDOS (not shown 
here) is negligible on the molecule, which explains their 
similar zero-bias conductance. 

Finally, the self-energy is also modeled by the MPM [^, 
[HI, which includes off-diagonal elements, using A=2 eV. 
The MPM leads to an effective opening of the gap be- 
tween the unoccupied and occupied molecular-like or- 
bitals which is much bigger (~4 eV) than the one ob- 
tained with the CHSX approximation (~2 eV). As a re- 
sult, the MPM conductance is much lower than in any of 
the fully ah initio approaches. This discrepancy can be 
attributed to some limitations in both the fully ab initio 
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calculations and the model ones. On the one hand, the 
use of highly-packed geometries considerably increases 
the screening and reduces the size of the corrections ob- 
tained within the fully ab initio calculations. Never- 
theless, the use of the CHSX approximation, which is 
known to overestimate the experimental gaps, may par- 
tially compensates this reduction. On the other hand, 
two main drawbacks can be pointed out for the model cal- 
culations: (i) the classical image charge model is not al- 
ways appropriate in predicting the individual shifts [2^; 
(ii) many-body effects on the gold region as observed in 
our calculations within the CHSX approximation. 

In summary, we have demonstrated that a self-energy 
operator leading to the mere correction of the eigenen- 
ergies (such as the GqWq or CHSX approximations) is 
not enough to change the initial DFT zero-bias conduc- 
tance. In contrast, when updating also the wavefunctions 
(such as obtained by the CHSX approximation), the con- 
ductance is reduced, improving the agreement with the 
experiment. The conductance change can be attributed 
to both a reduction of the molecular character and an 
increase of the 63(^^2) gold character. Finally, the dis- 
agreement between the MPM and ab initio calculations 
originates from the value of the scissor-operator shift A 
and from the wavefunction changes on the gold atoms. 
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